Efimov physics from the functional renormaUzation group 

Stefan Floerchinger • Sergej Moroz • Richard Schmidt 



Abstract Few-body physics related to the Efimov effect is discussed using the functional renormalization 
group method. After a short review of renormalization in its modern formulation we apply this formalism to 
the description of scattering and bound states in few-body systems of identical bosons and distinguishable 
fermions with two and three components. The Efimov eff'ect leads to a limit cycle in the renormalization 
group flow. Recently measured three-body loss rates in an ultracold Fermi gas of ^Li atoms are explained 
within this framework. We also discuss briefly the relation to the many-body physics of the BCS-BEC 
crossover for two-component fermions and the formation of a trion phase for the case of three species. 
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1 Introduction 

The physics of two or three non-relativistic particles that can scatter of each other or form bound states can 
be well described by quantum mechanics. Nevertheless, we present here an approach to these problems that 
is based on a method borrowed from the modem formulation of quantum field theory. Why should this be 
useful? There are several reasons. First, to better understand few-body physics it is desirable to have several 
different approaches at hand. Intuitive pictures and quantitative results can be compared and improved. 
Second, by applying powerful but also involved methods to rather simple problems one can also learn 
how to improve these methods. Also, by establishing a unified description of different physical phenomena 
(possibly on very different energy scales) one can hope for fruitful interference between different fields. 
Some universal aspects might even be transferred directly. Finally, for the field of ultracold quantum gases 
one has the very concrete advantage that functional renormalization can not only be applied to few-body 
problems but allows in a very similar way also the description of many-body phenomena. 

Our formalism uses the quantum effective action r[tl/] which is a functional of the field i//{t, x) repre- 
senting in our case non-relativistic particles. Formally speaking, the quantum effective action /"[(/'] is the 
generating functional of one-particle irreducible Feynman diagrams. This means that tree level expressions 
for scattering amplitudes are exact when the vertices and the propagators are obtained from functional 
derivatives of /"[(/']. In some sense the effective action r[iff] is similar to the microscopic action S [iff] but 
with the difference that all quantum and statistical fluctuations are included. Variation of r[if/] with respect 
to tfr yields exact equations. The functional form of r[i//] and the connection to scattering and bound state 
properties will become more clear below. Here we remark that for a generic interacting field theory the 
calculation of /"[i^] is non-trivial. A method for this task is provided by functional renormalization. 
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A particularly interesting feature of few-body physics relevant for experiments with ultracold quantum 
gases is universality U. Broadly speaking this means that for a large class of microscopic properties (such 
as the interaction potential between two particles) the macroscopic, low-energy properties are independent 
of the microscopic details and depend only on a few parameters. In particular if the scattering length is large 
compared to the effective range of the interaction potential the macroscopic properties can be equivalently 
calculated from a contact potential (with a vanishing range). In a renormalization group picture universality 
manifests itself as a (partial) fixed point. A large class of microscopic theories is attracted towards the same 
RG trajectory during the flow towards the infrared. The low-energy properties of the whole class can be 
characterized in terms of one representative which is for simplicity chosen to be the contact potential. For 
this case the microscopic or classical action has a particularly simple form. For a single boson field ij/ (as 
appropriate to describe identical bosons) it reads 



The microscopic interaction strength will later be related to the s-wave scattering length a. We use in 
Eq. ([TJ a formulation with imaginary time r which is integrated from to 1 / T in the Matsubara formalism. 
We also introduced a chemical potential /y. Although we are interested here only in the few-body properties 
in vacuum, i.e. for vanishing temperature T - Q and particle density n - 0, the generalization to many- 
body physics at nonvanishing temperature is possible and straightforward in principle by adjusting fi and 
T. Real-time properties can be obtained from analytic continuation t - it. Throughout this paper we will 
use natural units with Ti - ks - 1 . In addition we measure energies in units of momentum squared which 
leads to 2M^ - 1 . 

Signatures of universality also arise in few-body physics related to the Efimov effect f^. A remarkable 
feature which underlies the Efimov effect is the presence of a discrete scaling symmetry. Specifically, 
at the unitarity point, where the scattering length diverges, the solution of the Efimov problem remains 
invariant (self-similar) under rescaling of length by integer powers of some discrete scaling factor In the 
renormalization group language the discrete scaling symmetry is reflected in a renormalization group flow 
that periodically traverses a closed orbit, a so-called limit cycle Q. Limit cycles encode the essence of 
the Efimov effect and broaden the concept of few-body universality introduced in the previous paragraph. 
Recent experiments with ultracold atoms (for a summary see iQ) provided the first clear evidence of how 
Efimov physics is realized in nature. 

Since the three-body sector develops a renormalization group cycle scaling at the unitarity point, one 
must provide a proper ultraviolet initial condition for the three-body interaction, which is dictated by the 
details of the microscopic physics such as for instance the underlying interatomic potentials. Phrased in 
RG language, the limit cycle scaling implies that a three-body term ~ Ai(iff*if/)^ has a complex anomalous 
dimension and is exactly marginal at unitarity. For these reasons a three-body interaction term must in 
principle be included in the microscopic action ([!}. In the main part of this work we do, however, not 
concentrate on the direct comparison with particular experiments in which case the microscopical details 
and correspondingly a three-body term would become of relevance. For the sake of simplicity we tune the 
three-body parameter to zero which is why such a term is not included in the microscopic action ([TJ. In 
this simple setting all three-body correlations are generated during renormalization group evolution. Only 
in section|6] where we compare our findings with recent experiments with three-component fermionic ^Li 
atoms, the three-body term becomes a part the microscopic description. 

This paper is organized as follows. In section |2] we give a short introduction to the functional renor- 
malization method with emphasis on the non-relativistic few -body problem. In section [3] we discuss the 
problem of two identical bosons with contact interaction as well as the problem of three identical bosons 
at a broad Feshbach resonance where the scattering length diverges. The subsequent section|4]is dedicated 
to a system of fermions with two species (e.g. hyperfine states). We discuss how the renormalization group 
flow equations are modified compared to the case of bosons in the few-body sector and sketch how they 
account for the BCS-BEC crossover physics at non-zero density. Section |5] contains a discussion of our 
renormalization group treatment for a system with three fermion species. We assume that the fermions 
have equal properties (apart from their hyperfine spin) and that there is a common Feshbach resonance for 
all particles leading to a global SU(3) symmetry. We generalize this model to the case where the resonances 
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occur at different values of the magnetic field in section |6] In this setup we explain three-body loss rates 
recently measured in experiments with ^Li atoms ||5]|6l. Finally, we draw our conclusions in section Q 



2 Method 

2.1 Functional renormalization group 

Functional renormalization is a modem realization of Wilson's renormalization group concept which is 
used in quantum and statistical field theory, especially when dealing with strongly interacting systems (for 
reviews see l?)). The method combines functional methods of quantum field theory with the physically 
intuitive renormalization group idea. The main motivation for its development was the observation that 
it is often more useful and transparent to perform the integration of quantum fluctuations in continuous 
steps rather than doing it at once. The renormalization technique allows to interpolate smoothly between 
the known microscopic laws and the complicated macroscopic phenomena. In this sense, it bridges the 
transition from simplicity of microphysics to complexity of macrophysics. 

In quantum field theory a particular useful quantity is the effective action F which is a direct analogue of 
the classical action functional S . It depends on the fields of a given theory and includes all quantum fluctua- 
tions. For a generic interacting field theory the effective action F, however, is difficult to obtain. Functional 
renormalization provides a practical tool to calculate F employing the concept of renormalization group. 

The central idea of functional renormalization is the introduction of a scale-dependent effective action 
functional Fk often called average action or flowing action. The average action Fi^ interpolates smoothly as 
a function of the RG sliding scale k between the known microscopic action S and the full quantum effective 
action F. The dependence on the sliding scale k is introduced by adding an infrared regulator term to the 
microscopic action. This regulator term has the purpose to suppress fluctuations below the scale k and it is 
introduced by defining (for a bosonic field i^) the ^-dependent Schwinger functional WkU} 



The difference to the standard QFT Schwinger functional W[J] is the additional term AS it[if/] in the expo- 
nential in Eq. AS is quadratic in the fields and in momentum space it reads 

ASkW^ f <P*(Q)Rk(Q)m). (3) 

Jq 

We use the abbreviations Q = (^o,q) and = with = ^ dqo and { ^ jd'q. The 

infrared regulator function Rk{Q) is not fixed in the unique way, but must satisfy three important conditions. 
First RkiQ) should suppress infrared modes below the scale k, i.e. 

Rm * k^. (4) 

This adds effectively an additional masslike term to the particle suppressing its propagation if its momentum 
is below the RG scale k. Second, the regulator must vanish in the infrared 

^mRk(Q) = 0. (5) 

With this condition it is guaranteed that Wt equals the standard, exact Schwinger functional W in the 
infrared, W - Wa=o- Finally, we should recover the classical action in the ultraviolet and thus we demand 

hm RkiQ) = oo. (6) 

Apart from these requirements, the regulator can be chosen arbitrarily. Nevertheless, it is recommended to 
pick the regulator carefully depending on a concrete physical problem. Moreover, one should better choose 
a regulator that respects as many symmetries of the studied problem as possible. For our purpose it is 
convenient to choose RkiQ) - Rkiod independent of the frequency argument ^o- Often, Rki(() is chosen to 
decay for large » As we have seen, in Eq. (|2]) the regulator Rk suppresses fluctuation modes with 
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Fig. 1 Renomialization group flow in the theory space of afl possible couplings allowed by symmetries. Different trajectories coiTe- 
spond to diff'erent choices of the infrared regulator Ri^. 



momenta < by giving them a large "mass" or "gap", while high momentum modes are not (or only 
mildly) affected. Thus, Wa includes all fluctuations with momenta > k^. 

The flowing action is now obtained by subtracting from the Legendre transform of the Schwinger 
functional, 

hW = j Jiff-WdJ] (7) 

with (A = 6Wk[J]/6J, the cutoff" term 

rkW^hm-AS^m. (8) 

The reason why the flowing action is the central object to study is that it obeys the elegant and 
exact functional flow equation 

dknW = ^STrdkRk (rf V] + RkT', (9) 

derived by C. Wetterich in 1993 [jS]. Despite its compact form, the flow equation governs the full quantum 
dynamics. In Eq. (|9]l dk denotes a derivative with respect to the sliding scale k at fixed values of the fields. 
The functional differential equation for ^[i/'] must be supplemented with the initial condition A^a - 
S, where the "classical action" S[i//] describes the physics at the microscopic ultraviolet scale k = A. 
Importantly, in the infrared limit ^ — > the full quantum effective action r[iff] - Fk^oliff] is obtained. In 
the Wetterich equation STr denotes a supertrace operation which sums over momenta, frequencies, internal 
indices, and different fields (taking bosons with a plus and fermions with a minus sign). 

It is important to note that the exact flow equation for Fk [tfr] has a one-loop structure. This is a significant 
simplification compared for example to perturbation theory, where all multi-loop Feynman diagrams must 
be included. Note however, that Eq. (|9]i is a nonlinear functional differential equation. Indeed, on the right- 

(2) 

hand side in the denominator appears the second functional derivative T^. [i/r]. Due to that it is usually not 
possible to find closed solutions to Eq. ^ in an analytical form. 

To make the implications of ^ more transparent we write the flowing action as an expansion in some 
complete set of operators 

oo 

n=0 

The renormalization group evolution of A traces a trajectory in the theory space, which is the infinite - 
dimensional space of all possible couplings {c„) allowed by the symmetries of the problem. As schemat- 
ically shown in Fig. \T\ at the microscopic ultraviolet scale k - A one starts with the initial condition 
fk^A = 5 . As the sliding scale k is lowered, the flowing action Fk evolves in theory space according to 
the functional flow equation (|9]i. Accordingly, different couplings {c„} become running or flowing under 
the renormalization group evolution. The choice of the regulator Rk is not unique, which introduces some 
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scheme dependence into the renormalization group flow. For this reason, different choices of the regulator 
Rk correspond to the different paths in Fig. [1] At the infrared scale A: = 0, however, the full effective action 
Fk^Q - r is recovered for every choice of the cutoff and all trajectories meet at the same point in the 
theory space. 



2.2 Non-relativistic few-body problem 

As was mentioned above we use here a field-theoretic formalism with which also many-body problems 
in thermal and chemical equilibrium can be addressed. In this setup the properties of a few particles such 
as scattering amplitudes and binding energies are encoded in the n-point correlation functions in vacuum. 
The projection of the effective action F onto the vacuum state was developed in ll9l [T0l . It amounts in 
cooling the system (T — > 0) and making it dilute (« — > 0) in such a way, that the temperature T always 
stays above its critical value in order to avoid many-body effects (e.g., Bose-Einstein condensation). The 
chemical potentials of different particles in a given problem are related to their binding energies measured 
with respect to some reference energy level. As usual in non-relativistic physics the absolute energy scale 
can be shifted, and it is often convenient for us to fix it such that the lowest energy excitation (lowest energy 
n-particle bound state) has a vanishing energy. 

All few-body systems discussed here are translation invariant, i.e. no trapping potential is considered. 
Moreover, they are invariant under space-time Galilei transformations, which restricts the form of the non- 
relativistic propagator of a field if/ to he a function of the combination dr - of temporal and spatial 
derivatives (we recall that in our units 2M^ - 1). In addition, all systems possess an internal U(l) symmetry 
associated with particle number conservation. 

As was discussed in detail in fToHTTHTSHTSl . numerous mathematical simplifications appear in few- 
body problems compared with more complicated many -body problems. For instance, all Feynman diagrams 
with loop lines pointing in the same direction necessarily vanish. This simplification originates from the 
fact that there are only particles, but no antiparticles or holes as excitations of the non-relativistic vacuum. 
As a result, a remarkable hierarchy of RG flow equations arises in few-body physics. The flow equations of 
the n-body sector are not influenced by the higher-body sectorsQ One can solve the two-particle problem 
without any information about the three-particle problem and so on. We note that the presence of this 
hierarchy makes the systematic vertex expansion, discussed in Sec. 12.3.2] very useful for the treatment of 
few-body problems. The generally applicable strategy thus is to solve n-body problems subsequently by 
starting with the lower (one-particle, two-particle) sectors and proceeding to the higher ones. 



2.3 Truncations 

In most cases of interest the flow equation (|9]l can only be solved approximately. Usually some type of 
expansion of F^ is first performed, which is then truncated at finite order leading to a finite system of 
ordinary differential equations. The choice of a suitable expansion scheme should be physically motivated 
and depends on the given problem. The expansions do not necessarily involve a small parameter (like an 
interaction coupling constant or 1 /N) and thus they are, in general, of nonperturbative nature. Here we 
describe two different expansion schemes that we used for studying few-body quantum problems with 
functional renormalization. 

2.3.1 Derivative expansion 

A particularly simple truncation scheme is the derivative expansion. One expands the flowing action in 
terms with increasing number of spatial and temporal derivatives^ As an example consider bosons with the 

' The n-body sector is defined here as the set of one-particle irreducible 2«-point functions expressed in ternis of the elementary 
i^-fields. For a detailed discussion of the few-body hierarchy as well as a general proof of this feature see 1131 . 

^ Due to the non-relativistic nature of the problem we can assign to the derivatives the scaling dimensions [V] = 1 and [dr] = 2 
which corresponds to the dynamical critical exponent z = 2. 
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microscopic action given in Eq. ([T]i. In the derivative expansion the ansatz for the flowing action is 

r [uk{ip*>P) + Ziip* drip + Z2ip*{-v^)iij 

Jt,x^ (11) 



Similar to the effective potential Uit{4f*4r), being a function of t//*ip, the coefficients Z, and V, are functions 
of the scale parameter k and depend in principle also on the U(l) invariant combination t//*if/. Besides higher 
order derivatives also terms of the form i/r*i/r(-V^)"i/r'i/r will appear in Eq. ( fTTT i. 

The derivative expansion is particularly useful for the investigation of many -body problems jT]. In 
problems where bound state formation plays a role one expects that the derivative expansion in the form 
(fTTT i is not particularly useful. The reason is that one expects in this case a pole in the effective four-point 
correlation function. Extracting this object from Eq. (ITTl by taking the appropriate functional derivatives, 
one finds that it is given by a combination of terms with different frequency and momentum dependence. 
However, all these terms are regular and any finite number of them cannot account for a bound state pole. 

For this reason, instead of working with the microscopic action in Eq. ^ it is sometimes useful to 
perform a Hubbard-Stratonovich transformation. One adds a term that is quadratic in the complex auxiliary 
field d). The action becomes 



Jt,x ^ 



(12) 



Since the action for (p is Gaussian, the field <p can be easily integrated out. Hence, the action (fTZt is equiv- 
alent to the original action ([TJ, and under the condition 

+ 0, (13) 

^0 

the interaction term ~ is cancelled. The resulting theory is of the Yukawa type 

5 = J {-A*(5, - - |/)-A - hifipip + -AW) + fmltp^- (14) 

It is possible to investigate the theory in this form using a derivative expansion for the original field tp as 
well as the composite auxiliary (dimer) field <p. This will be described in more detail below. 

For some purposes it is useful to do the Hubbard-Stratonovich transformation not at a single scale but 
in a A:-dependent way. This is advantageous in particular when the interaction parameter that develops a 
bound-state pole is due to fluctuation effects and emerges from the functional renormalization flow. An 
exact flow equation that realizes such a scale-dependent Hubbard-Stratonovich transformation has been 
derived II14II15I . For a discussion of this approach to the treatment of bound states in the context of non- 
relativistic physics see ref. |[T6l . 



2.3.2 Vertex expansion 

As we advocated in Sec. 12.21 a systematic vertex expansion of is especially suitable for studies of few- 
body problems. The vertex expansion is an expansion in powers of fields, i.e. schematically 

oo 

A = ^ A(«) = A(2) + A(3) + A(4) + (15) 

«=2 

where n represents the number of fields in the monomial Fkin). We dropped A(0) in ( fTSl l which encodes 
the ground state energy of the non-relativistic (Fock) vacuum. The term linear in fields is absent due to 
U(l) symmetry. It is important to stress that, even if absent in the microscopic action S - rk=A, the vertices 
/"/;(«) become in general nonlocal in position space during the RG flow. Assuming translation invariance it 
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is convenient to Fourier transform to momentum space, where the monomials A(n) are functions of « - 1 
kinematically independent momenta and energies. At the infrared scale {k - 0) this kinematic dependence 
of directly translates into energy and momentum dependence of scattering and bound state properties of 
the system. 

To make this construction transparent we consider the example of bosons described by the microscopic 
action S in Eq. ([1]). The term A(3) vanishes in this case, and the vertex expansion in Eq. (fTSl ) starts with 
the terms 

If . . 

A(4) = - MQi'Q2^ QiWiQiWiQiMQimQi + Q2- 

^ JQ1.-Q4 

where P^{Q) is the scale-dependent inverse propagator and A^{Qi, Q2, Q3) is the nonlocal scale-dependent 
four-particle vertex. Due to the U(l) particle number symmetry only terms with equal number of i/f and iff* 
can appear in the vertex expansion. 



3 Identical bosons 

In this section we investigate the few-body physics of identical bosons interacting via a contact interaction. 
The hierarchy of flow equations introduced in section IZ21 turns out to be extremely useful as it allows us 
to tackle few-body problems step by step. Completely detached from higher n-body sectors we can first 
analyze the flow equations describing a single particle, then proceed to the two-body sector followed by 
the three-body system. For solving the «-body problem we only need the solution of the m-body systems 
with m < n. 



3.1 Two-body problem 

The interacting bosons are governed by the microscopic action in Eq. ([TJ. In the following it will be useful 
to switch to the formulation after the Hubbard-Stratonovich transformation in Eq. (fT4l i. 

In order to apply the functional renormalization group equation (|9j we first have to choose a truncation 
for the effective flowing action F^- We choose a vertex expansion as introduced in Eq. (flST l. which, up to 
third order in fields, reads in momentum space 

r,(2)= r ^\Q){iqQ + q^ -^i^)l|f{Q)+ r 0*(0^0(e)<^(e), 
Jq Jq 

A(3) = ^ r [^'{Qi)m2)mi) + mi)riQi)r{Q3mQi-Qi-Qi)- (i7) 

We stress that A(2) and A(3) do not have the most general form compatible with Galilean and global U(l) 
particle number symmetry. This would include a general inverse atom propagator P^{Q) and a momentum 
dependent Yukawa coupling h{Q\, Q2). However, using the special properties of the vacuum one can show 
that neither P^{Q) nor hiQuQi) get renormalized. Indeed, the condition of vanishing density implies 
that the poles of the propagators with respect to the frequency argument ^0 are restricted to the upper 
complex half -plane. Therefore diagrams with all internal lines pointing in the same direction (we call these 
diagrams 'cyclic') have integrands whose poles lie in the upper complex frequency half-plane. By closing 
the integration contour in the lower half-plane one finds that the corresponding expressions vanish. 

Let us investigate now the renormalization of the two-body sector Here, we have to consider the cou- 
plings h, A^, and the inverse dimer propagator P^. At the UV scale, k = A, the contact coupling vanishes 
because it has been replaced by the dimer exchange interaction in the Hubbard-Stratonovich transforma- 
tion. Yet there is the box diagram corresponding to the y6-function d,A^, schematically shown in Fig.|2](left). 
This diagram potentially regenerates A^ during the RG flow. It is however cyclic and thus one has A^ = 
during the whole RG flow. The only possible diagram renormalizing the Yukawa coupling h contains the 
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dt 




Fig. 2 The beta functions of the couphng (left) and the dimer propagator (right) depicted schematically as one-loop Feynman 
diagrams. The arrows in Feynman diagrams correspond to flow of particle number associated with U(l) symmetry. 



vertex A^, which is zero. For this reason also h is not renormalized, and one is left with the flow equation 
for the inverse dimer propagator P^{Q). For this reason, the truncation (ITTt is complete to third order in 
fields. 

The diagram renormalizing is shown in Fig. |2] (right) and can be written as followsjl 



with p - I for bosons i/rQ 



3 +p J 



(P^(L) + R^{L)){P^{Q -L)+ R^iQ - L)) 



P^Q) = iqo + q - 



(18) 



(19) 



and R^iQ) denotes the atom regulator In Eq. (fTSl l. the symbol d, represents a formal t derivative that hits 
only the cutoff R^{Q). For the choice of R^ we follow IIIOIII 111 and choose a frequency and momentum 
independent cutoff of the simple form R^ = k^. It turns out that this masslike regulator allows a relatively 
simple, analytical integration of the flow equation (fTST l. The two-body problem has also been solved ana- 
lytically using the so called optimized Litim cutoff by Birse [|17|| . In contrast to the regulator, the Litim 
regulator breaks Galilean invariance such that symmetry restoring counterterms have to be inserted in the 
UV. 

First we perform the frequency loop integration in Eq. (ITSl with the help of the residue theorem: 



3+p J (2;r)3 iq^ + 



2h^ ^ f d^l 1 
377 'J (2n)^iqo+r^ + {l-q, 



q)2 - 2p^ + 2R^ 
q)2 - 2fi^ + 2R^ 



(20) 



In the second line we replaced di with dt- The tilde indicates that the derivative acts only upon the k- 
dependence of the regulator We are allowed to drop the tilde as neither the Yukawa coupling h nor the atom 
propagator is renormalized such that only the regulator term exhibits a A:-dependence. This permits to 
perform the t integration of the flow equation from the UV scale fuv = to the IR scale tm = -oo. Using 
the specific values of the regulator R^ihR) - and /?^(fuv) - ^ we obtain 



Pf{Q)-P''^{Q) 
2h^ 



2h- f d^l ( 1 
~J (2^1/?o + l2 + (l-q)2- 

/j2 r di 

~~3+p j 2^ 



1 



2^l^ iqo + f + q)2 - + 2A'- 



, . [iqo . 



(21) 



^ Here we use the logarithmic, dimensionless RG scale ; = In ^ such that di = kdk- 

* The scattering physics of two distinguishable particles such as two fermions with opposite spin can be solved analogously with 
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In the first line we get rid of the angular integration by substituting 1 — > 1 - q/2 and the last identity assumes 
A » |ju^|,|q|,M- 

At this point we have to fix the initial condition for P^^ at A: = yl in order to obtain the physical inverse 
propagator P^^ at A: = 0. The momentum and frequency independent part of the inverse dimer propagator 
at the microscopic scale can be written as 

=(5v + v(fl). (22) 

The physical properties determined by the inverse infrared propagator P^^, should be independent of the 
cutofl" scale A. From Eq. (1211 1 we see, that this condition determines the term 6v, 

= irr; r^- (23) 

An{2> +p) 

Furthermore, as the composite dimer mediates the interaction between atoms, P^^ has to correctly de- 
scribe the scattering physics characterized by the s-wave, zero-range scattering amplitude 

fiq) = (24) 

-a ^ - iq 

where E - 2q~ is the total energy of the two scattered atoms with momenta +q respectively. The scattering 
amplitude is in turn connected to Pjj^. A simple tree level analysis shows 

1 

m = -« = TSTTT ^ ' (25) 

Ibn P]l^(0,0,^^ = 0) 

where the dimer propagator is evaluated for vanishing frequency and momentum. At this point it is im- 
portant to note, that when describing the s-wave scattering of atoms one deals with atoms at the energy 
threshold in the in- and outgoing states. For this reason it is necessary to work at vanishing atom energy, 
fi^ = 0. Taking this into account we are in the position to determine v(a). Evaluating Eq. (1211 1 for vanishing 
momentum, energy and atom gap fj.^ we find 

v(a)^-— -a-\ (26) 

47r(3 + p) 

We finally arrive at the solution of the two-body problem encoded in the exact infrared inverse dimer 
propagator 

This propagator also recovers - after performing the Wick rotation to real energies - the exact (s-wave) 
zero-range scattering amplitude in Eq. (l24t . 

Quantum mechanics predicts the existence of a universal, weakly bound dimer state for two particles 
with an interaction characterized by a large, positive scattering length a. Its binding energy can be extracted 
from a pole of the T-matrix at negative energy. In our approach the T-matrix is given by 

t{E) = (28) 

P;>,0,/i^ = £/2) 

and finding its pole therefore corresponds to solving ^'^'^(0, 0, /j.^) - for the atom gap -jj.^ > 0. This leads 
to 

a = ——. (29) 
As the dimer consists of two atoms we find for its binding energy 

2 

Eu = 2/1^ = — Eu = -. (30) 

The second equation is expressed in conventional units and is the well-known universal relation for the 
binding energy of the shallow dimer The determination of //^j, takes care of the correct vacuum limit. The 
atom gap yu^;, < is fixed to the value where the dimer becomes gapless in the infrared. This reflects that 
for fl > the dimer is the two-body ground state of the system. 
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Fig. 3 Kinematics of the vertex /isCQ'^, Q'^,E) in the center-of-mass frame. The atoms and dimers have momenta g'f = (E^i,<ll) . 
et = (-E^il + E, -qi) and = (E^l.qi), Q\ = (-£^2 + E, -qz). 



3.2 Three-body problem 

In order to treat the three-body system of identical bosons within the vertex expansion scheme, we extend 
our truncation (fTTI i by the term 



A(4) = - j M{QuQi,Qi)mmQ2)4>\QiW{QA)6{Qi + Q2-Qi-Q,) (31) 

21,-24 



with At,{Q\, Q2, Qi) denoting the atom-dimer vertex. For simplicity we restrict our discussion to unitarity, 
fl ' = 0, and demonstrate how Efimov physics emerges from the renormalization group flow equation for 
the atom-dimer coupling A^iQx, Q2, Q^). 

First consider the kinematics of the atom-dimer scattering. The IPI vertex At,{Q\, Q2, Qi) depends on 
three independent energies and spatial momenta. In the following we work in the center-of-mass frame and 
take the incoming atom and dimer to have momenta qi and -qi, and energies E^i and E - E^\, while the 
outgoing particles have momenta q2 and -q2 and energies £^2 and E-E^2- In the following the atom-dimer 
vertex will be denoted by /i3(2p 62' The kinematics is summarized in Fig. [3] 

The closed exact solution of the two-body problem, presented in the previous subsection, provides a 
simple strategy for the derivation of the renormalization group flow equation for the atom-dimer vertex. 
In principle, one can introduce completely independent regulators and for the atoms ijj and dimers 
(f>. Nevertheless, since the presence of the cutoff" R^ does not alter the solution of the two-body problem, 
we choose for the atoms the Galilei-invariant momentum independent cutoff" R^, used in subsection 13. II 
and lower it first from R^ = to zero while keeping R^ fixed. As the result of this step the dimer inverse 
propagator is renormalized according to Eq. (l27i . This also induces interactions between dimers, as for 
example a term ~ {(p*(p)^. However, these interactions belong to the four-body and higher sectors. By virtue 
of the vacuum hierarchy, they do not influence the flow of At,. Subsequently in the second step, only the 
dimer cutoff' is present and we lower R^ from its UV value to zero. In the second step and are fixed 
according to Eqs. (IZTT i and ( fT9l l. In this subsection we use a sharp cutoff' for the dimer field. 



(32) 

The special feature of this cutoff is that the regularized dimer propagator takes the simple form 



^ =«(lq|-^)ir77^- (33) 
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Fig. 4 Graphical representation of the flow equation for the atom-dimer vertex A3 . Full lines denote atoms ifr and dashed lines dimers 
0. The shaded circle denotes the atom-dimer vertex A3. 



Following the procedure described above we derive the flow equation for the atom-dimer vertex Aj(Q'^, Q^; E), 
which is given by 



P^{L)P4-L + Q) 



^1 

^2! 



CA,{QlUE)A,(L,Qt;E) 



+A- 



A^{L,qI-E) + A,{Q\,L-E) 



(34) 



P^{-L+Q\) P^(-L+Ql) 

where Q - + Q\ - (E, 0), A = 1, B = 2 and C - I. The graphical representation of this flow equation 
in terms of one-loop Feynman diagrams can be found in Fig.|4] 

Fortunately, at low energies and momenta the flow equation can be considerably simplified (for details 
see mi). After introducing the angle-averaged IPI atom-dimer vertex 



Aiiquqi 



1 r' 



(i(cos0)/l3(qi,q2;£'), cos(0) = 



qi q2 

Iqillqal 



(35) 



and transforming to the real-time formalism, the flow equation reads 

2(3 + p) 



d,A-i(,q\,q2,E) 



CAj,{quk, E)A-i(k, q2,E) + 



I {^3(^1 , k, E)G{k, qo) + G(qu k)A3{k, qi, E)} + AG{qi , k)G(k, q2) 



where we remind that for bosons p - 1, and the symmetric function G(q 1 , q2) is defined by 

1 



Giquqi) 



4qiq2 



log 



q] + ql + qiq2 - j - ie 



q\ + q\ - qiq2- j - ie 



(36) 



(37) 



The infinitesimally positive ie term arises from the inverse Wick rotation and makes both Eqs. ( [36] l and 
(iJTl i well-defined. Note that Eq. (l36l l is independent of the Yukawa coupling h and thus is well-defined in 
the broad resonance limit /? — > 00. 

The flow equation ( |36] |. obtained above, is quite complicated. In order to gain physical understanding, 
we first employ a crude, but simple and intuitive pointlike approximation. In this approximation the IPI 
vertex Ai{qi,q2, E) is replaced by a momentum and energy independent coupling /I3. In the pointlike limit 
the flow equation ( l36b takes the simple form 



d,AS 



4(3 + p) 



(38) 



where the rescaled dimensionless coupling A^ - A^k^ is introduced. The right-hand side of this diff'erential 
equation is a quadratic polynomial with constant coefficients. The sign of its discriminant D determines 
the behavior of the solution of Eq. (l38T l. For the bosonic problem D x -7.762 and we obtain a periodic 
solution of the form A^ ~ tan[;rf/r] with a period T - -^=- Note that the renormalization group evolution 
is governed by a limit cycle scaling. RG limit cycles are the manifestation of the Efimov effect. Indeed, we 
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can give the following intuitive interpretation of the periodic solution; During the RG flow we encounter 
three-body atom-dimer bound states, which manifest themselves as divergences of the atom-dimer coupling 
/l^. In the unitary limit there are infinitely many of these bound states, which are equidistant in a logarithmic 
RG energy scale. In this way the continuous scaling symmetry is broken to the discrete scaling subgroup, 
and the generalized universality emerges. 

It is well-known that at the unitarity point, where the scattering length a diverges, the energy spectrum 
of Efimov trimers forms a geometric series 

^ = exp(-27r/^o) (39) 

with £,,+1 and E„ denoting neighboring bound state energies. The Efimov parameter io is given by the 
solution of a transcendental equation and one finds sq ~ 1.00624 |[T|. By dimensional arguments we can 
connect the artificial RG sliding scale with the binding energy E as E ~ k'. The proportionality factor 
disappears in the ratio of the energies and hence the Efimov parameter can be read off from the RG period 

-j^ = = exp(-2r) 50 = -. (40) 

In the pointlike approximation we obtain sq ~ 1.393, which diifers from the correct result by 40%. Ob- 
viously the simple pointlike approximation is too crude to get the correct quantitative agreement. Never- 
theless it provides us with the qualitative picture of how the Efimov eff'ect appears as limit cycles in the 
functional renormalization group framework. 

It is remarkable that for E = the flow equation ( l36l l can be solved exactly. Here only the main results 
will be summarized while for details we refer to fllOillll . In order to find the exact solution most easily we 
perform the redefinition 

f,(tut2,E) = 4{3+p)qiq2A3{quq2,E) gihJi) = 4(3 + p)qiq2Giquq2) (41) 

and introduce the connected, one-particle reducible atom-dimer vertex 

Mh , f2) = pftih , t2) + gih , t2), (42) 

where from now on we prefer to work with logarithms of momenta f i = ln{q[/A) and f2 = In{q2/A). The RG 
scale dependence of the atom-dimer vertex ft{t\,t2, E) is denoted by the subscript t. The renormalization 
group flow equation and the initial condition for /, are given in matrix notation by 

dj, = --^f, ■ A, ■ f„ ^0 = 8, (43) 

where A, has matrix elements A,{ti,t2) - S{t - t\)6{t - t2) and matrix multiplication denotes integration 
over /. The formal solution of this differential matrix equation reads for re (-oo, 0) 

where / denotes the identity matrix. 

In the IR limit t — > -oo, which corresponds to the integration of all quantum fluctuations, / = fi~-oo 
solves the following matrix equation 

f^g+-^g-f. (45) 

This is the well-known integral equation for the amputated, connected Green's function obtained by Skor- 
niakov and Ter-Martirosian more than fifty years ago [ITSl . 

For illustration purposes we discretized Eq. (l44t and solved it numerically. A series of cartoons of 
the RG evolution of the IPI vertex f,{t\ , t2) is depicted in Fig. |5] In the UV we start with a simple initial 
condition fi{t\,t2) - and observe how a periodic structure (with period Tspatiai ~ 6.2) gradually develops. 
This intricate periodic landscape is completely missed in the pointUke approximation which corresponds 
to a planar landscape with no f i and t2 dependence. The RG evolution of the point f,(ti = 0, fa = 0) of the 
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Fig. 5 The RG evolution of the momentum dependent modified vertex /f(?i,?2) = 4(3 + p)qiq2A^(q\,q2.E). Spatial momenta t\, 
t2 and the RG time t are discretized to = 50 intervals with a step At = 0.4. The cartoons correspond to the descritized steps 
10, 25,35,50. 



100 



50 ■ 




-50 ; . 



-100' 

-8 -6 -4 -2 

t 

Fig. 6 Numerical evolution of the point f,(ti = 0, fi = 0) with the RG time t. Renormalization group develops a limit cycle with a 
period a 3.1. 

landscape is depicted in Fig. |6l Similar to the solution obtained in the pointlike approximation we observe 
"temporal" oscillations. The period of these oscillations Ttemp ~ 3.1 which leads to the Efimov parameter 
So = * 1 .0. This is in a good agreement with the "exact" result io ~ 1 -00624 obtained from the 

^ Icmp 

quantum mechanical treatment and is only limited by our numerical accuracy. 



4 Two-component fermions 

4. 1 Few -body problem 

The solution of the two- and three-body problem, presented in the previous section, can be straightfor- 
wardly extended to a system consisting of two fermion species with equal masses M^^ = M^^_ near a broad 
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Feshbach resonance described by the Euclidean action 

5 = J - V2 - + fp^<p - /i(0V,iA2 + r2'^\(l>)]. (46) 

Here, the two species of elementary fermionic atoms if/i , i//2 are described by Grassmann-valued fields, and 
is a composite bosonic dimer with the inverse propagator P^. Wherever equal indices appear summation 
is implied. Note that in Eq. (|46] | the dimer (p is just an auxiliary field which has no own dynamics, i.e. its 
inverse propagator = P(i,.k=A is a constant. This fermionic system has a global SU(2) x U( 1 ) internal sym- 
metry with {iff I , i//2) transforming as a doublet and as a singlet of the SU(2) flavor subgroup. Two-species 
fermion systems ( |46] | near a broad Feshbach resonance have for example been realized experimentally 
using two hyperfine states of ^Li or '"'K ultracold atoms. 

The two-body problem can be solved analytically following the very same route that we have taken for 
bosons in subsection 13. II The only diff'erence for fermions consists in setting = -1 in Eqs. (fTSTl-dZTli. 
Similar to bosons, for a > there is a universal dimer bound state with the binding energy given by Eq. 

dsoli. 

The solution of the fermionic three-body problem is also directly analogous to the computation for 
bosons, presented in subsection 13.21 For two-component fermions one arrives at the flow equation (l34t 
with coefficients A = 1, B = -2 and C - I. The solution of this flow equation is, however, qualitatively 
very different from the bosonic case. This can be understood already in the pointlike limit, where the flow 
equation is given by ( l38l l. Crucially, for fermions the discriminant D is positive. We find D - 9.881, and the 

solution of the flow equation is of the form ~ tanh ^-^(f + /<•) j, where k is fixed by the initial condition. 

The solution in the pointlike approximation approaches fixed points in the IR and UV. The numerical 
solution of the full RG flow equation ( l34b . detailed account of which can be found in IfTTII . supports this 
finding. In summary, the renormalization group flow does not develop a limit cycle scaling and there is no 
Efimov effect in the three-body problem for two-component fermions with equal mass. 



4.2 Many-body physics 

Although the few -body physics of two-component fermions of equal masses is not as rich as for the case of 
identical bosons, this system exhibits interesting many -body physics 11911201 . We concentrate on the case 
of equal densities ^ — - . For small negative scattering length a ^ 0_ BCS theory becomes valid. 
Many-body quantum fluctuations lead to spontaneous U(l) symmetry breaking and a superfluid ground 
state. Also for small positive scattering length a — > 0+ the ground state will be superfluid, now given by a 
Bose-Einstein condensate of dimers. In between these two regions one finds a smooth BCS-BEC crossover, 
including in particular unitarity with a ' =0. 

The functional renormalization description of this system has been developed in the recent years |[T2l 
l2ni22l . It proved to be convenient to work with the Hubbard-Stratonovich transformed microscopic action 
in Eq. ( l46l l and to use a truncation of the flowing action A in terms of a derivative expansion. The most 
elaborate truncation studied so far Il22l reads 

-I- Uk{(p'(f>) - /!(0*iAiiA2 + <P4'l4^\) + A^^fcptfr'il/i^. (47) 

Here, the coefficients Zip, in^, Z^, A^, h, A^,/, and the efl'ective potential Uk{4'*4') ^re ^-dependent quantities. 

The truncation in Eq. (l47T i can be used both for the few-body physics (vacuum) and to describe thermal 
and chemical equilibrium. In the former case one finds the characteristic binding energy of the dimers 
as in Eq. ( l30t . Due to simplicity of the truncation it is also rather easy to calculate more complicated 
observables such as the scattering length of dimers qm- In the limit a — > 0+ the result qmIci - 0.59 Il22ll is 
in good agreement with the exact result om/o - 0.60 ||23]| obtained from solving the Schrodinger equation. 
The dimer-dimer scattering was also recently studied with functional renormalization by the Manchester 
group 112411251 . In ll25l all possible local coupUngs that can be generated in the fermion four-body problem 
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were included and it was found that aula - 0.58 + 0.2, where the error was estimated from the regulator 
dependence of the final result. 

At non-zero density the flow equations are modified, for example the particle-particle diagram renor- 
malizing the boson propagator in Fig. |2]has a stronger eff'ect. As a result for small temperatures the effective 
potential Uk{(p*(p) becomes of a mexican hat form at the macroscopic scale k -Q. This leads to spontaneous 
U(l) symmetry breaking and superfluidity. For a more detailed discussion of the flow equations as well as 
for the results of this approach we refer to ref . ni2JI22jE6l . 

5 Three-component fermions with global SU(3) symmetry 

As was discussed in section |4] there is no Efimov effect for two fermion species with equal masses. The 
main reason for this is Pauli blocking. An obvious question is what happens in the case of three species 
(flavors) of fermions. The question is of particular interest as systems with three hyperfine states have been 
realized in recent experimants |5j|6j. In this section we will consider the case where all fermion masses are 
equal as well as the three pairwise scattering lengths such that the system possesses a global SU(3) internal 
symmetry. The generalization to the experimentally accessible system of ^Li atoms, where the pairwise 
scattering lengths are not equal, is discussed in section |6] 

5.1 Exact solution at unitarity: Vertex expansion 

For the case of SU(3) symmetric fermions we extend the action (|46] | to 

5 = J {iA; [dr - - ^) + f. [dr - v2/2 + m^) 0,- + '^eijk {<p]>pj^k - (t>irjr^ ]■ (48) 

Eq. ( l48l l is the analog of a two-channel model where the composite bosons 0, have a non-trivial propagator 
already in the microscopic action. The analog of the single channel model in Eq. (l46T l corresponds to the 
limit of broad resonances — > oo, wxj — > oo. The inverse propagator of the bosons is then dominated by 
the constant part m^. For finite h the action in Eq. ( l48l l can describe also Feshbach resonances of finite 
width. 

The three Grassmann valued fields i/r, for the fermions can be assembled into a vector if/ - , \l/2, i//}) 
and the bosonic fields (p, form a vector <p - {(pi,(p2,(p2) ~ (lAaiAs- 'A3'Ai> '/'i'A2)- The action possesses an 
SU(3) X U( 1) internal symmetry with if/ transforming as a triplet 3 and <p as an anti-triplet 3 under the SU(3) 
flavor subgroup. The three dimers mediate the pairwise scattering in the three possible s-wave scattering 
channels. 

What about Efimov physics in such a system? The answer to this question is in fact quite simple and 
can be understood with the following quantum mechanical argument: To identify the bound state spectrum 
one has to solve the three-body Schrodinger equation. In the case of fermions the bound state wave function 
has to be totally antisymmetric. The total wave function can be separated into the orbital and flavor part. 
If we take the flavor part to be antisymmetric, (e,j/t|/ > \ j > \k >), the orbital part must be symmetric and 
therefore has the same symmetry properties as in the bosonic case. As the orbital part solely determines 
the bound state spectrum in a quantum mechanical calculation, the bound state physics of SU(3) fermions 
is identical to the bosonic case leading to the same Efimov spectrum. 

A proper analysis of the renormalization group flow equations confirms these expectations. It is more 
complicated than for identical bosons since two different terms ~ (p* <pil/* tf/ similar to Eq. dSTT l are allowed 
by the symmetries. Specifically, the analog to Eqs. (fTTl i and (ISTT i reads 

We have suppressed here the nonlocal structure (momentum dependence) of the inverse propagator and 
the couplings A^a and AT,b for simplicity. Nevertheless, in the broad resonance limit — > oo, the solution 
of the flow equation at a ' = can be constructed in close analogy to the one presented in Sec. 13.21 In 
particular one finds also a limit cycle with the same Efimov parameter so ~ 1 .00624. For details we refer 
tolfin. 
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5.2 Approximate solution away from unitarity: Derivative expansion 



Although the foregoing discussion formally resolves the question about the Efimov effect for SU(3) fermions, 
we want to elaborate a bit more on this case. The solution of the flow equation in terms of the vertex ex- 
pansion works well (and is numerically exact) at the unitarity point with a ' - and £ = 0. However, one 
would like to understand also the region away from this point, especially since for experiments the physics 
directly at the unitarity point is hard to access. In this subsection we develop a simplified truncation of the 
flowing action in terms of a derivative expansion. 

Based on the ansatz ( l49l l our first goal is to find a 'minimal' truncation which still captures the essential 
features of Efimov physics. Due to fermion statistics the Efimov trimer is expected to be a SU(3) singlet. 
Only the term proportional to A^a in Eq. ( |49] l can have a bound state pole due to the combination of all 
three flavours. From this perspective, the term proportional to Ay,, which plays a similar role as in the two 
component fermion case, is expected to be less relevant for the description of the three-body bound state. 
We expect it only to affect quantitative results but not the qualitative features and we will neglect it from 
here on. 

Instead of keeping the full momentum and energy dependence of and At,„ as in the vertex expansion 
(|49] l we will employ a derivative expansion 

A = r U-idr - V2 - + 4>* [a^ [dr - v2/2) + mj] cjy, + x* [A, [dr - V^/S) + m^^x 

Jt,x ^ 

^\^'ik {fi^j<Pk - (pi^'j^i) + 8 {fi^'iX - <Pi^iXl }■ (50) 

The inverse dimer propagator is now approximated by a gradient expansion in leading order in momentum 
and energy with a gap term and wave function renormalization A^. The truncation (ISOl l furthermore 
contains a single component fermionic field x which is a singlet under SU(3) and represents the totally 
antisymmetric combination i//ii//24'i- The Efimov trimer field which we also call sometimes the 'trion', 
mediates the original atom dimer interaction ~ A^a and is introduced by a Hubbard-Stratonovich transfor- 
mation analogously to the dimer field introduced in section 12.31 Similar to the dimer we approximate the 
inverse propagator for the field x in terms of a derivative expansion. 

The Yukawa-type term proportional to g describes how the composite field x is formed from its con- 
stituents. At the microscopic scale A; = A we choose g = 0, A^ = 1 and oo. The field is then an 
auxiliary field that decouples from the other fields. With this choice we recover the original microscopic 
action ( l48T l in the broad resonance limit. In the next subsection we will relate the initial values of the dimer 
gap mi{A) and the Yukawa coupling h to the scattering properties of the elementary fermions i^, . 



5.2.7 Two-body problem 

We proceed as in section|3]by subsequently solving the flow equations in the two- and three-body sector In 
the following analysis we will employ the regulator - r{k^ - p^)6{k^ - p^) with r = 1 for the fermions 
i/r,-, r = A^/2 for the bosons 0, and r = A^/3 for the trion field From the arguments presented in Sect. 13. II 
it follows that the fermion propagator and the Yukawa coupling h are not renormalized. The flow equation 
for the boson gap parameter reads 

and for the boson wave function renormalization we obtain 

Eqs. dSTT i. (|52] | can be solved by direct integration, for details we refer to ref. ||27| . 

The initial values for A^, m^, and h can be related to the scattering properties of the elementary fermions 
close to a Feshbach resonance. We can use the freedom to perform a rescaling of the field (pi and choose 
A^{A) - 1 . To fix the remaining parameters we note that the scattering length between the fermions is given 
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dtm^ ^ dt- >*-^ -{ ) dtg^dt 



Fig. 7 The beta functions of the three-body couplings and g depicted as Feynman diagrams. The dotted hnes with three arrows 
denote the composite trimer field x- 

by fl = where the couplings are evaluated at the macroscopic (IR) scale k - Q and for vanishing jjL. 

This fixes the initial value of to be 

mliA) = -—a-' + -^A - 2fi, (53) 

^ STT 07T- 

where the term -2yu accounts for the fact that the dimer consists of two fermions. 

The value of the Yukawa coupling h is related to the width of the Feshbach resonance. The scattering 
length close to a Feshbach resonance at magnetic field Bq and width AB can be parametrized by II 1911281 

a = ^ (54) 
B-Bo 

with ai,g the background scattering length. Furthermore, we may identify the boson gap mj,(k = 0, = 0) 
with the energy of the closed channel molecule, = i-Im(B - Bq), with jim being its magnetic moment (for 
more details see lITSl ). From this one can see that the Yukawa coupling li^ is proportional to the width of 
the Feshbach resonance: 

- -inabgHuAB . (55) 

Broad Feshbach resonances correspond to /z^ — > oo. However, we can study also the crossover from narrow 
to broad Feshbach resonances within the two-channel model Eq. (l48T l. 



5.2.2 Three-body problem and energy spectrum 
The flow equations of the three-body sector read 



2 6g^ A^k' 

d,m^ ^ — _ _ ^ ^ (56) 



7r2 (3A^^2 _ 2A^p + 2m2)2 



and 



2 2gh^ k^ {6A^k^ - SA^fi + 2mp 



9,g . m^d,a - i;^(,._^).(3^^,,_2A,^,2™2)r (57) 

While the expression in Eq. (|56] | and the second term in Eq. (|57] | are straightforwardly found from the 
diagrams shown in Fig.|7] the first term in Eq. (l57t needs further clarification. 

Although the original atom-dimer interaction ~ X^a has been replaced on the UV scale by the trion 
exchange using the Hubbard-Stratonovich transformation, an atom-dimer interaction term of the form 



I 



builds up again in the RG flow due to the box diagram shown in Fig. [8] As it turns out, it is essential to 
include this evolution in our RG treatment to recover Efimov physics. In order to do so, we use a scale 
dependent Hubbard-Stratonovich transformation, a method developed in 114111511 . By this we absorb the 
regenerated atom-dimer coupling X^^ into the flow of the Yukawa coupling g. 
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dfX 



Fig. 8 The box diagram leads to the regeneration of the atom-dimer veitex 
Stratonovich transformation it is absorbed in the flow of the Yukawa coupling j 



■l^i/, in the RG flow. Using a i-dependent Hubbard- 



As was shown in 014II15I . the Wetterich equation gets supplemented by an additional term from the 
scale dependence of the Hubbard-Stratonovich transformation. We choose this scale dependence such that 
the box diagram in Fig.|8]gets cancelled. This leads to 

/,4 ^5 {9A^k^ -lA^ii + 4ml) 
~ ~ lln'-g {k^ - fiy (3A^^2 - lA^fi + Imp^ ' ^^^^ 

This determines the first term in Eq. (l57t and with this choice A^^ s is enforced on all scales k. The scat- 
tering between fermions and bosons is now described exclusively by the exchange of the (scale dependent) 
trion bound state x- 

We now investigate the energies of the different states - elementary fermions i/r, ("atoms"), composite 
bosons (pi ("dimers") or composite fermions x ("trimers") - as a function of the model parameters (e.g. 
the fermion scattering length a). In general it is not possible to find a simple analytical solution to the 
flow equations ( l56b , ( fSTl ) and we have to resort on a numerical solution. To find the energy of the lowest 
excitation we adjust the chemical potential /i < such that the lowest excitation of the vacuum is gapless 
at ;t = 0. 

For a finite UV cutoff scale A we find that for a"' — » -oo, the lowest energy state is given by the atoms 
which implies fi - 0, > 0, and > 0. No bound states are present in this regime. 

In contrast, for a'^ oo, the dimer (p is the lowest energy state, // < 0, = 0, and > 0. In fact, we 
can use our analytical solution of the two-body sector to obtain the dimer energy. In the limit A/\fi\ oo 
we find 

^2 



= - 



32n V (327r)2 16;r 



(60) 



which gives the binding energy of the weakly bound dimer state Ed = 2jj.. For finite h and the choice 
A^{A) - 1 the dimer binding energy also depends on the "effective range" 

327r 

reff = (61) 

according to 

, 2 

2reff 



ED = 2fi^ — — 

'"eff 



a ^ 



(62) 



In the broad resonance limit (vanishing effective range) one recovers the well-known result Ed = -Ija^ - 
-n^KM^a^). 

Let us investigate to the regime of strong interaction close to the point with a"' = 0. In a region 
fl^j' < a"' < fl^2 , we find that the lowest energy is given by the trion, such that < 0, > 0, and - 
for A: = 0. In this region the lowest energy state is a bound state of three fundamental fermions ~ i//iij/2^3- 
The critical scattering length a^^ i<^c2^^ which the Efimov trimer becomes degenerate with the atom 
(atom-dimer) threshold, depends on the effective range. For narrow resonances /i^/A <sc 1 we find a linear 
increase with as a^} = -0.0013/i^ and a^} = 0.0085/;^. For broad Feshbach resonances, li^/A » 1, the 
depth of the trimer state becomes independent of the effective range and is determined solely by the UV 
cutoff scale a^^, a^,^ ~ A such that the cutoff scale can be viewed as the additional three-body parameter 
needed to determine the exact position of the universal Efimov trimer state H]. This crossover behavior 
from narrow to broad resonances is shown in Fig. |9](left). 
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In Fig. |9] (right) we show our result for the energy per atom //„ of the trimer as a function of at 
unitarity a"' =0. For narrow resonances we find fiu - -2.5 x lO^^/i"*. In the broad resonance limit h~ — > oo 
also jj.u depends only on the three-body parameter given by the cutoff scale, ju„ ~ A^. 

From the work of Efimov |'2l, and as we seen in section ISTI one expects a whole tower of infinitely 
many Efimov trimer states close to the resonance. In the following we will calculate this excitation spec- 
trum. In order to do so, we make the following observation. First consider the inverse propagators of the 
atoms, dimers and trimer after analytical continuation, r - it, to real frequencies 

P>l,/A^ - + Y - 2;u + v^ifi + W0/2 - p-/4), 

Px/'^x = -w^ + y - 3;u + v^ifi + - p-/9). (63) 

We note that and v^, can dependent only on the particular combinations of ju, a>, and p^ indicated above. 
This is a consequence of an additional symmetry of the real-time microscopic action S - under the 
time-dependent U(l) transformation i// — > e'^'tfr, (p — * e'^^'ip, and x — > e'^^'x if we also simultaneously 
change the chemical potential /i + £. In vacuum we do not expect any anomalies, such that also 
the quantum effective action F = possesses this symmetry. As consequence any dependence on // 
is accompanied by a corresponding frequency dependence and vice versa. This allows only the particular 
combinations co^ + jj., a>^ + 2fi, and + 3;U to appear. The additional relation between frequency and 
momenta is a consequence of Galilean invariance. 

The dispersion relations are obtained from the poles of the propagators, that is from solving P^(a)^ - 
S^) - and analogous for the atoms and the trion. The dimer dispersion for example reads 

= y - 2/" + v^(S^/2 + ^ - P'/4) = y + ^ ■ (64) 

To obtain the dimer and trimer energy levels and we have to substract two or three times the atom 
energy E^ = -fi from the dispersion relation at rest such that 

E4, = &4P^ = 0) + 2yu = v^cj^ = p^ = 0,fi), 

E^ = S^ip- = 0) -H 3;u = v^ioj^ = S^, p2 = 0, //). (65) 

Here ju is the effective chemical potential which accounts for an shift in the atom energy using the time- 
dependent U(l) symmetry. The ground and excited states are obtained by the condition that the particular 
particles are gapless in the IR which corresponds to a gapless dispersion, i. e. for example = 0) = 0. 

This means that the corresponding gap parameter has to vanish in the infrared, i. e. yu = 0, or nr - 0, or 
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Fig. 10 The Efimov spectrum for SU(3) fermions. We plot the lowest three Efimov trimer levels (blue). The dimer binding energy 
is shown as black curve. To improve the visibility of the energy levels, we rescale both the dimensionless inverse scattering length 
a^^ /A and the dimensionless energy E/A^. We have chosen the UV cutoff A = ag and the Yukawa coupling /rA = 100 as appropriate 
for 'Li atoms in the nif = 1/2, mf = -1/2 channel close to the Feshbach resonance at Bq ~ 810G. 



- 0. In practice, we extract the bound state energy levels - Ijj, = 3// by calculating the gaps 
and for k = such that they vanish for a specific choice of E^ - -ji. 

As it turns out for the trimer there is not only the ground state solution which we discussed before but a 
whole infinite set of possible solutions corresponding to the infinite tower of universal Efimov trimer states. 
The numerical result for the dimer energy level E^ and the lowest three Efimov trimer levels E^ as function 
of a"' measured in units of A is shown in Fig.[TO] In this figure we used a convenient rescaling to make the 
excited Efimov states visible. We have identified the UV cutoff" A with the inverse Bohr radius, A - ag\ 
and have additionally chosen h- - lOOA. This choice corresponds to the width of the Feshbach resonance 
of ^Li atoms in the nif - 1/2, nif - -1/2 channel relevant for experiments ||29l . 

5.2.3 Limit cycle scaling at the unitarity point 

Already in section|3]we have discussed the especially interesting unitarity point in the spectrum where the 
scattering length diverges and the atom energy is at the threshold, a"' - E^ - - Q. It is instructive to 
consider this limit also within the derivative expansion as it allows us to obtain analytical results. At the 
unitarity point all length scales drop out of the problem and one obtains a simple scaling solution for the 
flow equations. In this limit the flow equations for the three-body sector simphfy and read 

2 36 

^'"^ = 25^ ¥ 

where we used A^{k) - h^/{6n^k) and - A^k^ for k Q. Using g(A) = as initial coupling we find the 
oscillatory, analytical solution 

"Jt" = (j) iiiM)\<x>t\s« In ^) + sin (s„ In i J |, 

with So - V419/25 ~ 0.820 From this solution we see that tn^(A) drops out of the ratio g^/rn^ which 
determines the atom-dimer interaction A^a ~ g^/mi- The flow of and g^ describes a closed Umit cycle 



^ Using the analogous deiivative expansion for bosons, where there exists only a single term describing the atom-dimer interaction 
~ <p'i//'(/>il/, one finds .vq ~ 0.93 (30|. 
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with the frequency given by sq. Considering the simplicity of our approximation this is in reasonable 
agreement with the 'exact' result, io x 1.00624. The limit cycle goes on 'forever' only exactly at the 
unitarity point and at the threshold energy - 0. A finite a"' as well as a finite atom energy = -/i act 
as infrared cutoffs, which stop the limit cycle flow once k x \jj\ or k x 37r|fl"'|/4, respectively. 



5.3 Many -body physics 

The simple truncation of the exact functional renormalization group equation in terms of a derivative ex- 
pansion yields a detailed picture of the spectrum of bound states and scattering in the system with three 
species of fermions with SU(3) symmetry. Depending on the scattering length a, the lowest energy state 
may be a fundamental fermion (atom) or a bosonic molecule (dimer) or a composite fermion formed from 
three atoms (trion). Trions are stable (lowest energy state) in the interval (fl'i*)"' < a"' < ("3*2) 

We have so far only considered the flow equations in vacuum, but we can already infer some interesting 
features of what happens at non-zero density for the many-body ground state at T = 0. The qualitative 
properties for small density follow from our computations and simple arguments of continuity. Importantly, 
we find that a novel "trion phase" separates the BEC phase and the BCS phase 1271 . 

For nonvanishing density, the chemical potential /i increases above its vacuum value //o- At T = for 
fl"' > (flc2) ' the BEC phase occurs for small density. Due to the symmetries of the microscopic action, 
the eff'ective potential U{p) for the bosonic field depends only on the SU(3) x U(l) invariant combination 
p = A^(p^(p - A^{<p*(pi + (p*2(f>2 + iplipi)- It reads for small fi- fio 

t/(p)= yp'-2(yu-//o)p, (68) 

where we use the fact that the term linear in p vanishes for fi = fiQ, i.e. ni^ifio) - 0. The minimum of the 
potential shows a non-zero condensate 



with density 



Po = ^ , t/(po) = -— (ju-/Uo) , (69) 

dU(po) 4 

n = ^ = ~ J"o)- (70) 

For the BCS phase for a < Ud < 0, one has jUq = 0. Non-zero density corresponds to positive ju, and 
in this region the renormalization flow drives m^^ always to zero at some finite k^, with BCS spontaneous 
symmetry breaking (po > 0) induced by the flow for k < kc- Both the BEC and BCS phases are therefore 
characterized by superfluidity with a non-zero po = A^^ ' 0. 

As an additional feature compared to the BCS-BEC crossover for a Fermi gas with two components, 
the expectation value for the bosonic field ^ in the three component case also breaks the "spin" SU(3) 
symmetry. Due to the analogy with QCD this is called "Color Superfluidity" 1311 . see also ||32| . For any 
particular direction of the expection value of a continuous symmetry SU(2) x U(l) remains. According 
to the symmetry breaking SU(3) x U(l) — > SU(2) x U(l), the effective potential has five flat directions. 

For two component fermions the BEC and BCS phases are not separated, since in the vacuum either 
yuo = or »zj = 0. There is no phase transition, but rather a continuous crossover. For three component 

fermions, however, we find a new trion phase for (a^,j')"' < « ' < (flj.2) this region the vacuum has 
//o < and ni'^ > 0. The atom fluctuations are cut off by the negative chemical potential and do not drive 

to zero, such that for small density remains positive. Adding a term m^p to the effective potential 
(l68T l we see that the minimum remains at po = as long as m^^ > 2(/i - fio). No condensate of dimers 
occurs in the trion phase. The BEC and BCS phases, that show both extended superfluidity through a 
spontaneous breaking of the SU(3) x U(l) symmetry, are now separated by a phase where po = 0, such that 
the SU(3) X U(l) symmetry remains unbroken (or as we argue below will be only partially broken). 

Deep in the trion phase, e.g. for very small |fl '|, the atoms and dimers can be neglected at low density 
since they both have a gap. The thermodynamics at low density and temperature is determined by a single 
species of fermions, the trions. In our approximation it is simply given by a non-interacting Fermi gas. 
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with fermion mass 3M and chemical potential 3(iu - /jq). Beyond our approximation, we expect that trion 
interactions are induced by fluctuations. While local trion interactions ~ (x*X)^ forbidden by Fermi 
statistics, momentum dependent interactions are allowed. These may, however, be "irrelevant interactions" 
at low density, since also the relevant momenta are small such that momentum dependent interactions will 
be suppressed. Even if attractive interactions would induce a trion-trion condensate, this has atom number 
six and would therefore leave a Ze subgroup of the U(l) transformations unbroken, in contrast to the BEC 
and BCS phase where only Z2 remains. Furthermore, the trions are SU(3)-singlets such that the SU(3) 
symmetry remains unbroken in the trion phase. The different symmetry properties between the possible 
condensates guarantee true quantum phase transitions in the vicinity of a^i and for small density and at 
r = 0. We expect that this phase transition also extends to small non-zero temperature. 

While deep in the trion phase the only relevant scales are given by the density and temperature, and 
possibly the trion interaction, the situation gets more complex close to quantum phase transition points. 
For a Qci we have to deal with a system of trions and atoms, while for a a: 0^2 a system of trions and 
dimers becomes relevant. The physics of these phase transitions is not clear now, it might be complex and 
interesting. 



6 Three-component fermions without SU(3) symmetry 



Quite recently three-component Fermi gases have been realized with ultracold atoms l|5]|6l. In these exper- 
iments ^Li atoms are prepared in three hyperfine states. Using external magnetic fields the scattering length 
of the different pair of atoms (1,2), (1,3), and (2, 3) can be tuned due to the presence of three Feshbach 
resonances. In the SU(3) invariant system discussed before all pairs of atoms had equal scattering lengths. 
This is not valid in the case of ''Li atoms as can be seen from the scattering lengths profiles shown in Fig. 
nn Although the Feshbach resonances are close to each other, they are still at different magnetic fields and 
also the profiles of the scattering lengths differ for the different pairs of atoms. In order to compare our 
calculations with the experimental observations we have therefore to extend the model discussed in section 

m 

The experiments by Ottenstein et al. Q and Huckans et al. ^ concentrated on the magnetic field 
region between B - 100 G and B = 550 G where all scattering lengths are rather large and negative. The 
main observation was a distinct feature in the atom loss rate which could be attributed to a three-body 
process and shows up in the three-body coefficient which is defined by the loss rate: 

where n,- denotes the density of the various atom species. The measured data points from 13] are shown in 
Fig. [13] (right) together with our numerical results which will be derived below. 

Already from the work of Efimov f2\ it is known that trimer bound states should persist in systems 
were the pairwise scattering lengths of the three involved particles differ from each other as long as at 
least to of them are large in magnitude. This is the case in the regime investigated by the experiments and 
in the following we will show how the observed loss features can be explained by Efimov physics using 
functional renormalization 133] . Similar predictions have been made using different methods by Braaten et 
al. Il34l and Naidon und Ueda | 



6.1 Derivative expansion 

First consider the model which we will use to describe the system. We generalize the truncation given by 
Eq. (l50b presented in section|5] The average action then reads 




i.j,k=l i=l 
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Fig. 11 Two-body scattering lengths ai2 (solid), aij (dashed), and oij (dotted) in dependence of the magnetic field B of *Li atoms in 
the lowest three hyperfine states. These curves were calculated by P. S. Julienne 1291 and are taken from Ref. (5). 



The generalization consists in allowing for three different inverse boson propagators determined by the 
six couplings A^, and in^., i - 1, 2, 3, and the three independent Yukawa couplings /i,. This enables us to 
describe the distinct scattering properties of the three channels. Furthermore we allow for three couplings gi 
which couple a dimer and a fermion to the trimer bound state As in section|5] the trimer field mediates 



the effective interaction between dimers and atoms Afj such that for vanishing center of mass momentum 
one has 



^ = (73) 

By permitting three independent couplings g, we are able to account for different strength of the atom- 
dimer interaction. Note, however, that we neglect the atom-dimer scattering channel denoted by AT,b in Eq. 
(|49] l. We again expect its effect to be subdominant. 



6.2 Flow equations and initial conditions 

The flow equations of the two-body sector are the straightforward generalization of Eqns. (ISTT i and (l52t with 
the replacements — > — > A^,, and h — » h,. Since the Yukawa couplings /j, are not renormalized 
the flow equations can be immediately integrated. As in section|5] the microscopic values of the gap terms 

and the hj are chosen to reproduce the experimentally measured scattering lengths fli2,fl23,fli3 in the 
infrared. Due to our ignorance of the microscopic details (e. g. effective range) leading to the specific 
scattering profiles in the regime between B - 100 G and B - 550 G (see Fig. fTTT l. we choose here large and 
equal values of h - hi - !i2 = /13. This corresponds to contact interactions at the UV scale. However, we use 
the dependence on the specific choice of li to estimate the error within our truncation. We emphasize that in 
this problem the bosonic fields <pi are not associated with the Feshbach molecules of the nearby resonance. 
They should rather be seen as true auxiliary Hubbard-Stratonovich fields representing the T-matrix of the 
atom-atom scattering. 

The flow equations of the three-body sector read for the gap term, 

= f -0 , , 7 (74) 

and for the Yukawa-like coupling gi one finds 

g2h2hik' (ek^ - 5/i + ^) g3hhik' (ek'- - 5;u + ^) 

dtgi = ^ 2 — r~2- ^'^^^ 

3n^A^2{k^ - fi)^ (3^2 - 2// + ^) 37T'-A^3(k'- - (3^2 - 2// + 

The flow equations for ^2 and ^3 can be derived from Eq. (iTST i by permuting the corresponding indices. 
There is an additional contribution to the flow of gi arising from the box diagram which regenerates the 
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atom-dimer interaction during the RG flow. As explained in the previous section it can be computed using 
a ^-dependent Hubbard-Stratonovich transformation and the flows of the gi receive the additional contribu- 
tion 



.2 



Here we define (/, 7, /) - (1, 2, 3) and permutations thereof and find 



(76) 



d,A^ = (77) 

67r2A^,(^2-//)W-2;U+^)2 

Although this contribution is essential to obtain the full tower of Efimov states at the resonance it plays only 
a subdominant role for the properties of the lowest Efimov state. For simplicity we will therefore neglect 
this term in the computation of the loss coefficient ^^3. 

In the broad resonance limit, h, — > 00, Efimov physics is universal in the sense that it is only dependent 
on the given scattering lengths 0,^ and an additional three-body parameter UJ. In section |5] we found that 
this parameter may be given by the cutoff scale A. Here, we however use the ratio of the couplings g, 
and at fixed UV scale A as the three-body parameter, which is equivalent to using A as the parameter 
itself. Also we choose g = g\ - gi - gi at UV scale such that our three-body parameter is given by 
^(3) = -g^{A)/in^{A). Note, however, that the couplings gi evolve differently during the flow towards the 
IR. 



6.3 Numerical results and comparison with experiment 

We solve the flow equations ( l74b and dTSl l numerically. We note at this point that the zero-range limit is 
valid if the magnitude of the scattering length is much larger than the range of interaction given typically 
by the van der Waals length lydw, which for ^Li is lydw ~ 62.5 aq. Although the loss resonances appear for 
|fl,y| > 2ly^iw such that the zero-range approximation might be questionable, we will assume its validity. We 
will, however, estimate the error by using various values for the Yukawa coupling h, which is related to the 
effective range, as we have seen in section |5] 

The trimer energy level can be calculated in the same way as in the case of SU(3) symmetric fermions 
by determining the chemical potential fi such that = for k - Q. We find = for large enough 
values of 0,^ which indicates the presence of a trimer bound state. Now, the experimental observations offer 
a way to determine the three-body parameter inherent to Efimov physics: We simply adjust the initial value 
of /t*^*' such that the Efimov trimer becomes degenerate with the atom threshold exactly at the position of 
the observed loss feature at B = 125 G. This fixes the three-body parameter /l'^' to the physical value and 
we assume that it is constant for all other investigated magnetic field strengths B. 

Having determined the three-body parameter we are able to calculate the energy level of the trimer 
bound state as a first prediction. We find the existence of the trimer in the magnetic field range from B - 
125 G to B = 498 G. The binding energy £7 is shown in Fig.[T2]as solid line when choosing li^ = 100 Aq ' as 
appropriate for ''Li in the (1,2) channel. The shaded region is an error estimate obtained by choosing values 
h~ 6 (20 300 flg '). If we include the contribution due to the A:-dependence of the Hubbard-Stratonovich 
transformation given in Eq. (l76t we find a reduction of the trimer binding energy. The result is shown as 
dashed curve. 

As a second prediction we determine the loss coefficient K3 measured in experiments ||5]|6l. Within 
our model the observed three-body loss can be explained as follows below, an illustration is shown in Fig. 
[13] (left). So far we have shown the existence of the Efimov state in the regime of interest. This trimer 
state may the formed in a collision of three atoms. However, one must keep in mind that ''Li atoms have 
a complicated substructure with various bound states so that the composite trimer can become unstable. 
We account for this decay to not further specified degrees of freedom by a phenomenological decay width 
r^. For example, the decay might involve the molecular dimer states associated with the nearby Feshbach 
resonance. The trimer formation from three atoms proceeds via the exchange of the effective bosons 0,. 
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Magnetic Field B [G] 

Fig. 12 Upper panel: Scattering length ai2 (solid), 023 (dashed) and f/31 (dotted) as a function of the magnetic field B for *Li. Lower 
panel: Binding energy Ej of the three-body bound state ^=(4- 11^2 1^3- The solid line corresponds to the initial value /r = 100 c/q'. The 
shaded region gives an error estimate by using values of h in the range = 20 a^ ' (upper border) to = 300 ' (lower border). The 
dashed hne shows the reduction of the trimer energy when the refermionization of the atom-boson interaction is taken into account. 

Evaluating the matrix elements corresponding to the diagram shown in Fig. [T3] (left) we can estimate the 
loss coefficient to be proportional to 




The amplitude for forming a trimer out of three atoms is given by ^JLj higj/rn^.. The amplitude for the 
evolution from the initial state of three atoms to the unknown decay products further involves the inverse 
trimer propagator which in the real-time formalism reads 

We evaluate it in the limit of small momenta - (2, p,)^ — > 0, and small on-shell atom energies = p?, 
o) - 2, aii 0. That the incoming atoms in the recombination process are at the atom threshold also 
implies that we have an effective chemical potential jj - 0. The trimer gap in^ is then negative for B - 
125 . . .498 G. A further ingredient of the calculation of the matrix elements in Eq. ( fTSl l is the decay width 
which we first assume to be independent of the magnetic field B. Furthermore the evaluation of Kj, via 
Fermi's Golden rule involves unknown vertices and phase-space factors for the trimer decay such that our 
compuation contains an unknown multiplicative factor ck- In terms of p in Eq. (l78l l we obtain 

Ki^CKp. (80) 

The resulting curve is shown in Fig.[T3](right) and agrees reasonably well with the experimentally measured 
data points. To sum up, our estimation of Kj, involved three parameters: The location of the first loss 
resonance at Bq = 125 G fixes the three-body parameter. The overall amplitude of at the first resonance 
is adjusted by ck and the unkown decay width Fj^ is tuned to the width of the first resonance yielding a 
trimer lifetime of ~ 10 ns. The extension of the curve away from the first resonance then involves no further 
parameter. 

Our calculation predicts a second sharp resonance feature at B = 498 G where the trimer becomes again 
degenerate with the atom threshold. The assumption of a constant fixes the width of this peak. The 
measured feature is however much broader. This may have a simple explanation; It is plausible to assume 
that the trimer decays into the molecular dimer states associated with the nearby Feshbach resonances 
and a free atom. Approaching the Feshbach resonances the molecular state becomes closer to the trimer 
energy level. Therefore it is reasonable to assume that the matrix element for the decay from the trimer to 
the dimers becomes larger. This in turn implies a shortening of the lifetime of the trimer as the magnetic 
field B is increased. We have tested this assumption and indeed found a significant broadening or even 
disappearance of the second resonance peak. Recently, there had been phenomenological studies of this 
reasoning in 1361 . 
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Fig. 13 left: Thi'ee-body loss involving the Efimov trimer. right: Loss coefficient K3 in dependence on the magnetic field B as 
measured in (5] (dots). The sohd line is the fit of our model to the experimental curve. We use here a decay width that is 
independent of the magnetic field B. 



By now, the system of three-component ^Li atoms is quite well understood also away from the mag- 
netic field range discussed above. This includes theoretical predictions 07II38I as well as experimental 
measurements O9II40I which also finally led to the first direct association of a Efimov trimer in ultracold 
atoms by Lompe et al. RT\ . 



7 Conclusions 

During recent years non-relativistic few-body physics has experienced a revived interest and undergone 
considerable progress. A wealth of experiments with ultracold atomic gases together with a significant the- 
oretical effort initiated the (re)birth of an exciting field known as Efimov physics. In this paper we reviewed 
the studies of various aspects related to the Efimov effect using the method of functional renormalization. 

By constructing solutions to exact functional renormalization group flow equations, we have obtained 
a rather detailed picture of non-relativistic few-body problems related to the Efimov eff'ect. We considered 
identical bosons and fermions with two or three spin components, the latter both for the case of a common 
Feshbach resonance with a global SU(3) symmetry and for the case of ''Li where this symmetry is broken. 

Solving the flow equations in terms of a vertex expansion leads to (numerically) exact solutions. A 
remarkable hierarchy in the structure of the flow equations allows to construct the parts of the quantum 
effective action relevant for the one-, two- or three-body problem without further approximations. Beyond 
the present context it becomes apparent how the quantum mechanical formalism is embedded in non- 
perturbative quantum field theory. 

Besides the vertex expansion we also established the derivative expansion as a way to find approximate 
solutions. In connection with a A:-dependent Hubbard-Stratonovich transformation this allows for a simple 
and very efficient parameterization of the essential physics. No equivalent approximation scheme is known 
to us for the quantum mechanical formalism. As we have demonstrated this scheme is particularly valuable 
in more complex situations where a full solution in the vertex expansion becomes too involved. An example 
for this was given by the calculation of the three-body loss rate in a system of ultracold *'Li atoms which 
shows how our theoretical investigations can be directly tested in experiments. 

Both in the vertex expansion and in the derivative expansion the Efimov effect shows up as a limit cycle 
in the solution of the flow equations at the point where the scattering length diverges. This mechanism is 
very interesting from a theoretical perspective. Realistic models in relativistic quantum field theory that 
show a similar behavior have not been found so far Also the universality of the non-relativistic few-body 
physics is very interesting from a theoretical point of view. While in a renormalization group description 
universality is usually associated with the presence of fixed points, the Efimov effect with its limit cycle 
scahng and the connected discrete scaling symmetry goes beyond this standard paradigm. This generalized 
notion of few-body universality is expected to extend beyond the physics of three particles. First studies of 
four-body universality using functional renormalization have been undertaken I.30J . 

In conclusion, we beheve that the functional renormaUzation group is a powerful and suitable method 
for investigation of non-relativistic few-body physics and of the Efimov effect in particular It allows to 
address few-body problems from an unconventional perspective and provides an interesting alternative 
to the standard quantum mechanical and effective field theory treatments. Importantly, the method is not 
only limited to few-body problems but also is widely used for studying many-body physics. In this way 
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functional renormalization provides a unified framework for studying an interplay between few-particle 
and many-particle phenomena. 
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